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| We calculate the next-to-next-to-leading order correction to the cross section 

for top quark pair production in e + e~ annihilation in the threshold region, 
resumming all O [(a s /(3) n x (a 2 s , j3 2 , a s j3)\ terms of perturbation series. We 
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CO , find that the magnitude of the NNLO correction is comparable to the size of 
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the NLO corrections. 
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43 . I. INTRODUCTION 



The cross section of hadron production in e + e annihilation belongs to the best known 



quantities in high energy physics. Far away from quark thresholds the cross section is well 
approximated by the results obtained in perturbative QCD (for a review see [|IJ). The 
situation is not so clear at quark thresholds which, however, are known to be of importance 
for a number of physical applications. 

Among such applications, a special place is occupied by a threshold production of tt 
pairs at the Next Linear Collider. It was suggested in Ref. [Q], that the large width of 
the top quark provides a natural cutoff for long-distance effects and, therefore, reliable 
predictions for the tt threshold production cross section are possible in perturbative QCD. 
Since then, the threshold production cross section of tt was studied in great detail 0. 
The commonly accepted conclusion is that one can perform precision studies of various 
quantities of direct physical interest (top mass, top width, strong coupling constant, etc.), 
once accurate measurements in the threshold region are conducted. However, all these 
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studies were performed using predictions for the top threshold production cross section 
valid up to 0(a s ) and, therefore, suffered from the ignorance of higher order QCD effects. 

It is worth emphasizing, that calculation of radiative corrections to the threshold cross 
section differs from standard perturbative calculations, which are done for higher energies. 
The difference is because of the fact, that, close to the threshold, the conventional perturba- 
tion theory breaks down ||. The physical origin of this phenomena is known from quantum 
mechanics: considering Coulomb potential as a perturbation, one gets series in a/ (3, where 
a is the strength of the potential and j3 is the particle velocity. When the velocity is small, 
this ratio becomes large and meaningful predictions can be only achieved once the series is 
resummed. It was demonstrated in , that if such resummation is performed, the threshold 
cross section becomes proportional to the square of the Coulomb wave function at the origin. 
In Ref. this result was generalized to the situation, when the produced particles are un- 
stable. It was concluded, that the cross section for the quark pair production is proportional 
to the imaginary part of the non-relativistic Green function of the QQ system, evaluated 
for complex energies. 

Since then, it was also realized, that the 0(a s ) corrections can be easily incorporated, 
because contributions of soft and hard scales completely factorize with this accuracy. The 
absence of this factorization property, as well as the technical difficulties with explicit higher 
order calculations, were the stumbling blocks in achieving the 0(a 2 ) accuracy. It is remark- 
able, that new results, obtained in the last several years, permit a relatively easy determi- 
nation of these corrections. 

In what follows, we present the calculation of the threshold cross section for the ti pair 
production which is valid with 0(a 2 , a s (3, j3 2 ) accuracy. 

II. THE FRAMEWORK OF THE CALCULATION 

We first discuss a framework of our calculations and introduce all relevant notations. 
The threshold region is characterized by a small value of the quark velocity f3: 



/ 4m 2 . , 

/J = yi-_<1. (1) 

To order 0(a 2 , a s (3, /3 2 ), dynamics of slowly moving quark-ant iquark pair is governed by a 
non-relativistic Hamiltonian 1 : 

H = H + V 1 (r) + U(p,r), (2) 

B. = £-Z±, (3) 
m r 

Vi(r) = -^{2ftln(M'r)+a 1 

47IT 



1 One can describe the QQ system by means of the non-relativistic quantum mechanics to this 
order since the radiation of real gluons shows up only at C(/3 3 ) order. 
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In the above equations, the strong coupling constant is evaluated at the scale /x: 

a s = a s (n). 



(5) 



(6) 



The scale // equals to fie 1 , 7 is the Euler constant. 

The operator U(p, r) is the QCD generalization of the standard Breit potential [JTU |. The 
last term in Eq.([|) is the non-Abelian contribution, originating from a correction to the 
Coulomb gluon exchange, caused by a magnetic gluon [fTTjl . The potential Vi(r) represents 
a deviation of the static QCD potential from the Coulomb one. It was calculated to order 
12 and to order a 3 in ]y|]. The coefficients there read explicitly: 
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(7) 



For the SU(3) color group, the color factors are Ca = 3, Cp = 4/3, Tr = 1/2. N L = 5 is 
the number of quarks whose masses have been neglected. 

Given the Hamiltonian H, one can find the Green function for the Schrodinger equation: 



(H — E — iS)G(E; r, r') = 5 {3 \r - r'). 



(8) 



Once the Green function is found, the cross section of the non-relativistic QQ pair 
production in e + e~ annihilationP] is obtained as: 



o- Is) 



Ana' 
~3s~ 



R(s) 



(9) 



2 In what follows, we consider only photon mediated process and do not take into account the 
Z-boson exchange. The axial-vector coupling of the Z-boson contributes O{0 2 ) relative correction 
to the threshold cross section. The vector Ze + e~ coupling is also suppressed, but can be taken 
into account in the same way as the photon contribution, which we treat in this paper. 
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where 
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In Eq.(^), we have included the 0{f3 2 ) correction originating from the expansion of the 
vector current which produces and annihilates a heavy QQ pair in the triplet S"-state. The 
quantity R(s) will be the central object for further discussion. 

If a calculation of R will be attempted, one will find, that the Green function at the origin 
does not exist because there are terms in the Hamiltonian H, which behave as l/r n , n > 2, 
for small values of r. The difficulty originates from the fact, that the region r — > is not 
properly treated in the Hamiltonian. Indeed, small values of r correspond to a region in the 
momentum space, where a typical momentum transfer between Q and Q, is of the order of 
the quark masses and therefore quarks cannot be considered as non-relativistic. For this 
reason, the use of the Hamiltonian H in actual calculations leads to the divergencies, which 
appear for r — > 0. 

The way to circumvent this difficulty is as follows. In order to perform a calculation, 
one introduces a cutoff A, such that a s m C A <m. The momenta region where k <C A is a 
non-relativistic region and can be described using a Hamiltonian H. The momenta region 
with k ^> A is a relativistic one and the calculation in this region should be performed 
using the rules of quantum field theory. We note that this is rather standard procedure 
for calculations, related to bound state problems. It is also well known that its practical 
realization often requires substantial effort. 

However, there is a possibility to use the result of the non-relativistic calculation with the 
cutoff in the following way: one takes the limit of the obtained result, considering kinematic 
region where a s <C (3 <S 1. In this particular region, the non-relativistic results are still 
valid; on the other hand as long as a s //3 <C 1, the resummation of the Coulomb effects is not 
necessary. Therefore, in this particular region, one can calculate the corrections applying 
the standard rules of the quantum field theory. In the framework of QCD, such calculations 
have been performed recently in Ref. ||15|| . One therefore can match the result of the non- 



relativistic calculation with the cutoff, directly to the result presented in |TJj and in this 



way completely eliminate the cutoff dependence. This procedure was suggested in |TJ[] and 
we will call it a direct matching procedure, in accordance with that reference. 

In what follows, we pursue this program in QCD. We confine ourselves to a strictly 
perturbative approach and we do not attempt any discussion of non-perturbative effects. In 
order to accommodate the phenomenologically relevant case of the unstable top quark, we 
will consider the total energy E as the complex variable E — > E + iT t , in the spirit of Ref. 



We have to mention at this point that a consistent implementation of the effects related 
to the finite width of the quark is not attempted in this paper. To 0(a s ) order such effects 
were studied in M. A reliable treatment of these effects to 0(a 2 ) is not available at the 
moment. 
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III. CORRECTIONS TO THE GREEN FUNCTION AT THE ORIGIN 



Let us first consider the Hamiltonian H , Eq.(fJ), as the leading order Hamiltonian. The 
corresponding Green function will be denoted by Gc(f,r'). 

We first discuss a correction to the Green function Gc(r, r') at the origin, caused by the 
operator U(p,r) (see Eq.([5])). The first order correction is: 

- J d 3 r G c (r, 0) U(p, r) G c (r, 0). (11) 

As long as we are interested in the QQ pairs, produced in the triplet S-states, only a 
corresponding projection of the operator U (p, r) should be considered. Substituting S 2 = 2 
and L = into Eq.(H), we get: 

tti \ P 4 , ll7rC V a * x(3)/ \ C Fas fl 2 1 C A C F a 2 s 

At this stage, it is advantageous to express this operator in terms of the zeroth order 
Hamiltonian H in order to apply the equation of motion for the Green function Gc{v, 0): 

(H -E)G c (r,0) = 5^(r). (13) 

This is most easily done by using the commutation relations: 

rrr . . 47r5 (3) (r) 2L 2 C F a s 

[H ,ip r ] = ^ + — - -V 1 , 14 

11 2 4ri«M + _2^ 
r ) r m mr z 

where p r = —i(d r + 1/r) is the radial momentum operator. Using these relations, one finds 
that the operator U(p, r) from Eq.(|T2]) can be presented in the form: 

U(p r ) = _ 3CfCLs \ Hq I) + 11CFa ° [ Hq ip 1 _ (2^f + 3C A ) C F a 2 s 

Am Am I ' r J 12m ' r 6mr 2 



Let us consider the first three terms in Eq.(|16|). Inserting them into Eq. (pTTj) and using 
the equation of motion for the Green function (|13D, we find: 



p C F a 



H + 2ip r 

nn nr'< 



2m 2 mr 



G c (r',r'')+Jd 3 rGa(r',r){^+ 3 -^f-\G G (ry). (17) 



We note, that all terms in the above equation, which do not contribute to the imaginary 
part of the Green function have been omitted. 

Leaving aside the "surface" terms in Eq . (|17D , one sees that the discussed perturbation 
can be absorbed into the zeroth order equation flT3|). Therefore, the corresponding correction 
to Gc(r, r') can be taken into account to all orders by finding an exact Green function for 
the Schrodinger equation 
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(H-E)G(r,r') = 6®(r-r'), (18) 
with the modified Hamiltonian 



H=^ + V(r), V(r) = -^^1 + - + V l( r)-^-^ + ^), (19) 



and the modified eigenvalue 



E 2 

£ = E + —. (20) 
4m 



It is clear, that the solution of the above equation will deliver G(r, r'), which is definitely 
valid to NNLO accuracy. Moreover, such a solution provides a resummation of some second 
order corrections. 

But there is more important reason to believe that such a treatment is more appropriate 
for the subthreshold (-y/i < 2m) energy region, than the first order time-independent pertur- 
bation theory. Let us consider the energy region below the threshold. For stable quarks, one 
would observe an appearance of narrow resonances in this region. It is well known, that for 
the realistic value of the top quark width (r t ~ 1.5 GeV), the resonances are smeared. Still, 
the excitation curve exhibits a maximum close to the position of the lowest lying resonance. 

For the purpose of discussion, we write an expression for an exact Green function: 

where ij} n> k and E n ^ are exact eigenfunctions and eigenvalues, which correspond to a Hamil- 
tonian H . When we calculate this Green function using time-independent perturbation 
theory, we expand both numerators and denominators of the above equation in power se- 
ries. One readily realizes, that this procedure is not so harmless for energy denominators, 
especially when the energy E is close to the position of the resonance. On the other hand, 
if the Green function is obtained as the solution of the Schrodinger equation, one gets the 



result directly in the form of Eq. (pi]) . It is mainly for this reason, that we think it is more 
safe to solve the Schrodinger equation exactly, than to perform the first order perturbation 
theory calculations. 

It is not so straightforward, however, to obtain a numerical solution for such a problem, 
since the perturbation 1/r 2 (cf. Eq. (|T9"D ) is too singular at the origin. As was already 
indicated above, the proper approach is to introduce a cutoff ro and to extract all terms 
which have a non-analytic dependence with respect to ro- On the other hand, all terms 
which have extra powers of ro, so that the limit ro — > can be taken, will be set to zero. 
Later on, the non-analytic ro-dependent terms will be determined by a matching of the 



result of such calculation to its perturbative counterpart jyj. In the next section we will 
show, how this procedure can be implemented for the numerical solution of the Schrodinger 
equation (|TB|). 

The "surface" terms from Eq . (|i~7|) will be discussed later. We note here, that these terms 
are linear in Go(r, r') and, therefore, contain at most the first-order poles in energy. Hence 
they do not contribute to the shift of the energy levels and there is no need to account for 
them beyond the first order. 
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IV. NUMERICAL SOLUTION OF THE SCHRODINGER EQUATION 

In this section, we demonstrate, how numerical solution of the Schrodinger equation with 
the potential V(r) (cf. Eqs. (|l"8"D and (|19D) can be constructed. The Schrodinger equation 
for the .S-wave Green function is written as: 

' ; ' r 2( +V(r)-£]G(r,r') = ^—5(r-r'). (22) 



m 



h - — 

dr 2 r dr 



It is convenient to define a new function g(r,r r ) = rr'G(r,r'), so that the Schrodinger 
equation simplifies: 



— + m[£-V(r)}\g(r 1 r') 
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According to the standard rules, the Green function can be written as: 
g(r, r') = A [#<(?")#> (r')6»(r' - r) + (/)#>(?")% - r')] 



(23) 



(24) 



where <7<,>(r) are two independent solutions of the homogeneous Schrodinger equation which 
satisfy proper boundary conditions at r = and r = oo, respectively Also, the constant A 
is defined by the jump of the derivative of the Green function at r = r' . 

The solutions g<(r) and g>(r) are constructed using two other independent solutions 
of the Schrodinger equation, g±(r), with a prescribed behavior at the origin [|J. However, 
since the potential V(r) in the Schrodinger equation is as singular as 1/r 2 , one would have 
problems with setting the standard || boundary conditions for g± at r = 0. To overcome 
this difficulty, we extract the leading asymptotic of the functions g±(r) for small values of r 



9±(r) = {mC F a s r) d± f±{r) 



where 



- (l ± VT^Ik 
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(25) 



(26) 



Substituting g±{r) to the Schrodinger equation, one obtains an equation for the function 
f± (r) , which is now free of the 1/r 2 term: 
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m 



dr' 



+ 



2d± d 
r dr 



+ 



g + (i + Cl i n(;U V) 2 + c 2 ln(/i'r) + 

\ 



c 3 



/±(r) = 0. (27) 



The coefficients ci_3 can be easily obtained using Eq.(^) and Eq.(f|). One derives then the 
asymptotics of the function f±(r) for r — > 0: 
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1 - mC F a s rA(r) + 0(r 2 ) 
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A(r) = /ii ln 2 (^V) + h 2 ln(/i'r) + /i 3 , 
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(28) 
(29) 

(30) 

(31) 



In the above equations, / and d stand for f± and d±, respectively. 

There could be some doubts about the validity of the boundary conditions, as derived 
from above equations, since oL is of order a 2 s and therefore the second power of a s appears 
in the denominator. Without going into explanations at this point, we mention that upon 
careful inspection the above boundary conditions appear to be absolutely legitimate. Later, 
we will present more detailed arguments in favor of such conclusion. 

The Green function is constructed as follows. The solution g < (r) is identified with g+(r). 
The solution g>(r) is constructed from the solutions g± in such a way, that the boundary 
condition at the infinity is satisfied: 



9>(r) = g-(r) + B g + (r), g>(r)^0, 
Therefore, one finds: 



for r — > oo. 



B = — lim 

r— >oo 



9-(r) 



9+(r) 



(32) 



(33) 



If the potential V(r) is real, the imaginary part of the Green function is proportional 
to the imaginary part of the coefficient B. However, this is not the case for the present 
problem: the Coulomb part of the potential V(r) is energy dependent and we consider the 
energy to be a complex variable. Therefore, the formula for the Green function at the origin 
should be modified. The modification is however simple. It is obtained in a straightforward 
way from the available asymptotics of the functions ip±{r). The result reads 



Im G(r , r ) 



m 3 (Ci?a s 
AttW 



Im I (mCFa s r 



,2(24 



B — Aj 



(r ) - A-(r )} 



(34) 



where W = —mCFa s (d + — gL) is the Wronskian of two independent solutions of the 
Schrodinger equation. 

Let us comment on the role of the second term in the above equation (0). If the 
potential V(r) in the Schrodinger equation were real, the A±(r) would be real as well. 
Taking the imaginary part, one then completely removes the contribution of the second 
term in Eq.([34|). In our case E is a complex variable, and the coefficient h% has a non-zero 
imaginary part. Moreover, this imaginary part is formally of the order of T t /(mdJ), which 
should be considered as a contribution of order unity. 

However, the limit k — > should exist for the Green function constructed as above. 
Therefore, the purpose of the last term in Eq. (|34]) is to cancel the 1/oL singularities of the 
function B, obtained with the boundary conditions presented in Eqs.(p9|-|3~T|). We note in 
this respect, that, by switching off the logarithm-dependent perturbations in the expression 
for V(r), one obtains an exactly solvable Hamiltonian, so that the statements made above 
can be easily verified. The corresponding discussion can be found in Appendix. 

Eq. ( J53J ) then explicitly demonstrates, that the non-analytic dependence on the cutoff r 
is indeed extracted and all power corrections with respect to the cutoff are neglected. This 
precisely corresponds to the desired form of the Green function of the Schrodinger equation. 
The residual dependence on r is removed using the direct matching procedure as described 
in the next section. 
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V. MATCHING AND FINAL RESULT FOR R 



In Eg . (p4[) we have taken into account all corrections to the ratio R which are due 
to relativistic effects in the quark-antiquark interaction and which can thus be called the 
dynamic ones. However, we still have to consider kinematic corrections that are: i) 0(f3 2 ) 
correction to s; ii) 0(p 2 /m 2 ) correction to the quark current (cf. Eg. (JlQ|) ) and, finally, iii) 
the "surface" terms from Eg . (|17|) . 

Using eguation of motion for the Green function flT^D, we obtain the result for R at 
NNLO, with both types of corrections included^: 



Rnnlo = 7;N c e 2 Q C F a s [ 1 + C X C F ( -± ) + C 2 C F ( ^ 



^N c e 2 Q C F a s \\ + C x C F (^ 
x . = Im — — 

vr^ 7 ^ H 6 




K^fol^^B-^W-^W \. (35) 



Here we have factored out all energy-independent corrections. They are parametrized by 
the constants C% and C 2 , which are divergent in the limit r — > 0. For this reason we use 
ah = a s (m) as the expansion parameter for these "hard" corrections. 

To get rid of the r -dependence we use the direct matching procedure, suggested in Ref. 
TA\. For this we consider y/s > 2m, set the width of the top quark T t to zeroQ and equate 



our result fl3~5|) to its perturbative counterpart [|15| in the kinematic region a s <C /3 <C 1, 
where both are supposed to be valid. We also set \i = m, so that a s coincides with ah- 

Let us note, that the direct matching procedure fixes the linear combination of C 2 and 
ln(mro) 



a x 2 



7T 



C 2 C F [ — ) -2Kln(mr ). (36) 



If we were working strictly to NNLO, this last combination would be the only thing we 
need for the final result. However, because of the large difference in scales, which govern 
relativistic and non-relativistic physics, we would like to write Eq.(|35[) in a factorized form 
and include an exact dependence on ro into the non-relativistic Green function. For this 
reason, we have to set a factorization scale. We do this by choosing ro in such a way, that 
the correction to the Coulomb Green function due to the 1/r 2 perturbation in the region 
a s "C P 1 is given by log(/3), without additional constants (see Appendix for more details). 



3 We note that, strictly speaking, the "surface" terms in Eq. ( |l7| ) were derived for the Coulomb 
Green function, but we substitute the "exact" Green function instead of the Coulomb one in our 
final formulas. 

4 The strongest divergence in these factors is 1/ro, so it not obvious a priori, that the functional 



form of Eq.(34) can be preserved in Eq. fl35|) . Upon careful analysis, this turns out to be possible 



to NNLO, however. 
5 Note, that in this case the functions A±(ro) drop out from Eq.(|35|). 
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Any other choice of ro would correspond to other (also legitimate) value of the factorization 
scale. 

A factorized form fl3"5|) of our final result makes sense only if a dependence on a choice 
of the factorization scale is weak. We have checked that changing the value of the cutoff 
between vq/2 and 2ro for r$ given by Eq.(^), we obtain small (~ 1 — 3%) variation of the 
resulting values of R. 

Therefore, accoring to our choice, we fix the value of the cutoff 



and obtain finally: 
Ci = - 



Co 
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e 2 ^ 
2m 



CfC^ + CaC^ A + TrNlCo 
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where 
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a = — . 

2 9 



9 
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(39) 



Eq. (|35j ) with definitions provided by Eqs. ( p7| -|39|) is our final result for the top quark 
threshold cross section with the NNLO accuracy. 

For numerical purposes, we have chosen m = 175 GeV and T t = 1.43 GeV. As an input 
value for the strong coupling constant we used a s (Mz) = 0.118. Fig.l provides our final 
results for -Rnnlo as a function of ^/s — 2m in comparison with LO and NLO results, for 
three values of the soft scale \x = 50, 75, 100 GeV. One can see that the NNLO corrections 
are as large as the NLO ones. 

There is also a moderate scale dependence of the NNLO corrections in the vicinity of the 
resonance peak. The position of the resonance peak appears to be sensitive to the variations 
in the scale [i on the level ~ 100 MeV. We note in this respect, that the shift of the ground- 
state energy due to the Breit perturbation is well known (see Ref. [ID| and Appendix) and 
its expected variation with /x is close to this value. 



VI. CONCLUSIONS 

We have presented a calculation of the next-to-next-to-leading order corrections to 
the threshold cross section of the top quark pair production in QCD, summing all 
O [{a s/ (3) n x (a 2 , /3 2 , a s (3)] terms of the perturbation series. We have found, that the NNLO 
effects are quite sizable. 

We have also discussed how the numerical solution of the Schrodinger equation with a 
singular potential can be constructed. 
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Let us comment on the large size of the NNLO corrections. We have checked that taken 
separately, both the Breit perturbation and 0{a1) terms from Vi(r), provide comparable 
contributions of the same sign. When we take them into account simultaneously in the 
Schrodinger equation, the NNLO contribution gets enhanced by roughly a factor of two, in 
the vicinity of the resonance. 

When we were writing this paper, we received the preprint |16| , where the same problem 
was studied. Our qualitative conclusion about the size of the NNLO corrections agrees with 
the conclusion reached in ||16|| . 
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APPENDIX 

In this Appendix, we discuss the construction of the Green function G(£;r ,r ) in the 
exactly solvable model, which is described by the Hamiltonian 

H = P ^ + V(r), V(r) = -^(1 + c ) - " ( 40 ) 
where £ and c are complex parameters. 



The Green function is derived following our discussion in Sect. [TV]. It is convenient to 



introduce a new variable z = Cpa s mr for further discussion. The solution of Eq. (|27|) , which 



satisfies boundary conditions given by Eq.(28), can be written as 



f ± (z) = e^F (d ± - iu{l + c), 2d ± , -™ J , (41) 

where v = Cfcl s / (^y^/m) . Using the asymptotic form of the confluent hypergeometric 
function for Re x = Re(—iz/v) 1, 

F(a } b,x)^^\e x x a -\ (42) 
T(a) 

one obtains the coefficient B, using Eqs. fl25D and (j33|): 

u(j.\ _ 1 ( hA 2-2 d+ r(2rf_)r(^ + - iv{\ + c )) 

The functions A±(r) reduce now to the constants: 

AAr) = (44) 
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We therefore arrive at the final expression for the imaginary part of the Green function 
ImG(r ,r ) for this model (cf. Eq . (j34f) ) : 

Let us first demonstrate that the proper limit k — > exists for the imaginary part of the 
Green function defined through Eq.( [45|) . In this limit, our model reduces to the ordinary 
Coulomb problem, so that the Green function in Eq.(fl^) should give us the imaginary part 
of the Coulomb Green function at the origin. To see how this happens, we expand the Green 
function Eq.(H5|) in power series in k. The first term in the expansion of B(v) is equal to 
(1 + c)/(2k). This term is completely canceled by the last term in Eq. fl4*5|) . The next term 
in the expansion provides the ro-independent result: 

lim G(r ,r ) = — -lmH(u,p), 

where 

H{v, /3) = ^ - (1 + c) [ T + m(-i/3) + V (1 - HI + c))] , ijj{z) = A l n r(z), (46) 

which exactly coincides with the imaginary part of the Coulomb Green function at the origin. 

We then set c = and expand Eq.fl45|) up to O(k) to obtain the correction to the 
Coulomb Green function due to the 1/r 2 perturbation. The result can be written as: 

Im [5G(r , r )] = ^^lm {h(v, pf + -^-^ , (47) 

where the value of ro from Eq.(|37D has been used. For stable quarks, the last term in Eq. (f47|) 
can be disregarded. For unstable quarks, it contributes an 0(YQ/mQ) relative correction to 
the Green function at the origin, which is beyond the intended accuracy and can be omitted. 

Let us also emphasize one advantage of the Green function as obtained from Eq. (|45|) . 
Consider the stable quark case. Then, for negative energies, the Green function should 
deliver the first order poles which correspond to the appearance of the QQ bound states 
in the spectrum. Eq.([45"l) shows, that such poles are provided by the singularities of the 
function T(d + — ii>{l + c)). The corresponding eigenvalues of the Hamiltonian ( flD^ ) are 

m(C F a s ) 2 (l + c) 2 
Sn ~ 4(n-d_y ■ (48) 

Using the relations (cf. Eqs. ( [T9|) and (f20|)) 

£ = E + ^ c=|^, (49) 



4m 2m ' 



as well as 



d± 4(l ± VT^), K= <^(H + |), (50) 
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one easily finds that the energy levels are located at 



E 



m(CFa s ) 2 m{Cpa s ) 



An 2 



+ 



TV 



11 

64n 



1 

6 



C f 



4C f 



(51) 



We note, that this is precisely what one gets, if the energy shift due to the perturbation 
is calculated using the standard rules of quantum mechanics. 
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FIG. 1. i?LO (dotted lines), -Rnlo (dashed lines), -Rnnlo (solid lines) as a function of energy 
— 2m, GeV. In all three cases, three curves correspond to different choices of the soft scale 
= 50 GeV (upper curves), fj, = 75 GeV and \x = 100 GeV (lower curves). We also use 
= 175 GeV ; = 1.43 GeV and a s {Mz) = 0.118 as the input parameters. 
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